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Abstract. The 7-ray emission above 250 GeV from the BL Lac object Markarian 421 was observed by the CAT Cherenkov 
imaging telescope between December, 1996, and lune, 2000. In 1998, the source produced a series of small flares, making it the 
second extragalactic source detected by CAT. The time-averaged differential spectrum has been measured from 0.3 to 5 TeV, 

which is well fitted with a power law: ^ cc ^Tov • I" 2000, the source showed an unprecedented activity, 

with variability time-scales as short as one hour, as for instance observed during the night between 4 and 5 February. The 2000 
time-averaged spectrum measured is compatible with that of 1998, but some indication of a spectral curvature is found between 
0.3 and 5 TeV. The possibility of TeV spectral hardening during flares is also discussed, and the results are compared to those 
obtained on the other TeV BL Lac, Markarian 501. 
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1. Introduction 

Among Active Galactic Nuclei (AGNs), blazars are those 
radio-loud objects having a jet pointing towards the observer, 
which has a relatively high buUc Lorentz factor, giving rise 
to a strong Doppler boosting of the observed fluxes. Blazar 
emission is dominated by their jet power output, which is 
mainly non-thermal, extending over more than fifteen energy 
decades. At low energies, their featureless optical continuum 
as well as their strong radio and optical polarization are due 
to synchrotron radiation of the magnetized plasma jet. At 
high energies, blazar jets also show remarkable properties. 
During the last ten years, the EGRET detector, operating above 
30 McV on board the Compton Gamma-Ray Observatory, 
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definitively opened the field of high-energy astrophysics by 
revealing that most extragalactic strong 7-ray emitters were 
blazars (von Montigny et al. 1995; Hartman et al. 1999). Their 
7-radiation power often dominates their entire spectrum, and 
it must be produced in a small enough region to account for 
rapid variability, already observed on time-scales less than one 
hour at TeV energies ( [Gaidos et al. 1996 ). Although the origin 
of jets is still uncertain, the study of their 7-ray emission in 
blazars can shed light on the nature and content of their plasma 
(e+e^ pairs or e^p), and give new insight into the high-energy 
particle acceleration and cooling processes occuring at the 
sub-parsec scale. Ultimately this could give useful information 
on how jets take form and lead to a better understanding of 
the energy extraction mechanisms in the surroundings of the 
central supermassive black hole. 

Markarian 421 (Mkn 421) is the closest known BL Lac 
blazar (at a redshift of 0.031), and the first discovered in the 
high and very-high energy (VHE) ranges. It was first detected 
as a weak source by t he EGRET inst rument up to a few GeV 
during summer 1991 (Lin et al. 1992). Eight months later, the 
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Fig. 1. The Mkn 421 7-ray signal seen with the CAT imaging telescope, for the most significant observations (signal significance 
greater than 3 cr in ^^30 min) taken in 1997-98 and 1999-2000 at zenith angle (?z<45° (left panel), and for the highest flare, 
recorded between 4 and 5 February, 2000 (right panel). In each panel, the main plot shows the pointing angle, a, distribution for 
ON (solid line) and OFF data (points with error bars, same duration as ON data), and the bottom-left inset gives the "ON-OFF" 
distribution; events are selected by the Q4, Qtot, and 7'(x^) cuts (see text), and the values quoted in legend come from the 
additional cut q;<6° (first two bins). The bottom-right inset shows the significance map of event excess per bin of 0.05°, obtained 
from the "ON-OFF" distribution of reconstructed angular origins (no cut on a): the source lies at the centre, and the axes denote 
the relative right ascension and declination with respect to the source coordinates (i?As and (5s). 



Whipple Observatory detecte d a clear signal fr om this object 
between 0.5 and 1.5 TcV (iPunch et al. 19921). Since then. 



Mkn 421 has been confirmed many times as a VHE source 
by various atmospheric Cherenkov imaging telescopes, e.g. 
again by the Whipple Observatory (Krennrich et al. 1997), 
by the stereoscopic system of HEGRA (Petry et al. 1990 
and by the CAT (Cherenkov Array at Themis) experi- 
ment ([Piron et al. 1998|, |1999b[). Along with Markarian 501 



(Mkn 501), Mkn 421 has thus become one of the two extra- 
galactic sources of the Northern hemisphere which has been 
firmly established in the VHE range. It has been also one of the 
most studied blazars and the target of several multi-wavelength 
observation campaigns from the radio band to the 7-ray 



range (see, e.g.. 


Macomb et al. 1995 




Buckley et al. 1996 


Takahashi et al. 1996t 


rakahashi et al. 19991 


2000 


Chariot et al. 1998 Maraschi et al. 1999 


. Recently, a ma- 



jor step was achieved during the 1999-2000 winter by the 
CELESTE atmospheric Cherenkov sampling experiment, 
which detected Mkn 421 for the first time around 50 GeV, fil- 
ling the last energy gap still remaining on this source between 
satellites and ground-based detectors. These observations, 
made in part sim ultaneously with th e CAT telescope, have 
been presented by |^older et al. (200l")| . 



In this paper we concentrate on the temporal variabiUty and 
the spectral properties of Mkn 421 above 250 GeV, as seen by 
the CAT experiment since it began operation in autumn 1996. 
Sect. ^ describes the detector and the analysis methods used to 
extract the signal and the spectra. The Mkn 421 data sample, 
the light curves up to June 2000 and the corresponding spectra 
are presented in Sect. |[ We discuss our results in Sect. ^ com- 
paring them with those from other ground-based atmospheric 
Cherenkov telescopes, and with those obtained on Mkn 501. 
The conclusions are given in Sect. ^ 

2. Experimental setup and data analysis 

2. 1. Telescope characteristics 

Located on the site of the former solar plant at Themis in 
the French Pyrenees (2° East, 42° North, altitude 1650 m 



above sea level), the CAT telescope (Barrau et al. 1998) uses 
the Cherenkov imaging technique. Its 17.8 m'^ mirror collects 
the Cherenkov light emitted by the secondary particles pro- 
duced during the development of atmospheric showers and 
forms their image in its focal plane. The CAT camera has a full 
field of view of 4.8° but the present study is limited to the 3.1° 
fine-grained inner part, which is comprised of 546 phototubes 
with 0.12° angular diameter, arranged in an hexagonal matrix. 
This very high definition, combined with fast electronics whose 
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integration time (12 ns) matches the Cherenkov flash's short 
duration, is an efficient solution of the two main problems of 
night-sky and cosmic -ray backgrounds which confront such 7- 
ray detectors. Firstly, the detection threshold, which is deter- 
mined by the night-sky noise, is as low as 250 GeV at Zenitl{]. 
Secondly, the capability of rejecting the huge cosmic -ray back- 
ground (protons, nuclei) is improved by an accurate analysis 
based on the comparison of individual images with theoreti- 
cal mean 7-ray images. This method, very briefly reviewed in 
Sect. 2.2, is specific to the CAT experiment and has been des- 
cribed in detail in |^e Bohec et al. (1998^ . 



2.2. Gamma-ray signal extraction 

In order to improve the hadronic rejection and stabilize the 
background level near the detection threshold, and to compen- 
sate for possible slight changes in the detector response be- 
tween different epochs of observation, we eliminate the noi- 
siest pixels and require the fourth-brightest-pixel's charge in 
the image Q4>3 p.e. (photo-electrons) and the image's total 
charge (3tot>30p.e. 

An efficient discrimination between 7 and hadron-induced 
showers is then achieved by looking at the shape and the 
orientation of the images. Since 7-ray images are rather thin 
and ellipsoidal while hadronic images are more irregular, a first 
cut is applied which selects images with a "7-like" shape; it is 
based on a fit to a mean light distribution predicted from 
electromagnetic showers, and a probability 'P{x^)>0.35 is 
required. In addition, since 7-ray images are expected to point 
towards the source angular position in the focal plane whereas 
cosmic-ray directions are isotropic, a second cut a<6° is used 
in the case of a point-like source, where the pointing angle 
a is defined as the angle at the image barycentre between 
the actual source angular position and the source position as 
reconstructed by the fit. As a result, this procedure rejects 
99.5% of hadronic events while keeping 40% of 7-ray events; 
the Crab nebula, which is generally considered as the standard 
candle for VHE 7-ray astronomy, is detected at a 4.5 cr level in 
one hour 

Fig. |l] shows the a distributions obtained from two data 
samples taken on Mkn 421, for ON and OFF-source observa- 
tions (the latter being taken at the same telescope elevation in 
order to monitor the hadronic background), and the correspon- 
ding distributions for 7-rays obtained by "ON-OFF" subtrac- 
tion (bottom-left insets). The signal is clearly seen in the di- 
rection of the source (small a), though the direction of some 
7-rays is mis-identified, giving a small signal at Q!'^180°. As 
stated above, the fit also allows the angular origin of 7-ray 
events to be determined with good accuracy as it uses the in- 
formation contained in the images' asymmetrical longitudinal 
light profile. In Fig. |l], the bottom-right insets show the signifi- 



The threshold is defined here following the usual convention for 
atmospheric Cherenkov detectors, i.e. as the energy at which the dif- 
ferential 7-ray rate is maximum at the trigger level for a Crab-like 
spectrum. As shown in Fig. ^, it is somewhat higher after the event 
selection applied at the analysis level. 




E,(TeV) 

Fig. 2. CAT effective detection area (top) and differential trig- 
ger rate (bottom) for 7-ray showers after event selection. 
Top: each point in the main panel is a simulation, while full 
lines come from an analytical 2D-formula AcsiOz, E^) allow- 
ing interpolation over zenith angle and energy. The inset com- 
pares the prediction of this interpolation for 6*2=20° with an 
independent simulation (open circles); 

Bottom: differential 7-ray trigger rate (^ x ^eff ). for a typical 

'2.55, 



spectrum ^1=3.0 i^T^y x 10 
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angles (6lz>45°). 
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cance map of 7-ray event excesses: the angular resolution per 
event is 0.11° (i.e., of the order of the pixel size), allowing a 
bright source to be localized with an accuracy better than 1' 
(dominated by systematics). 

2.3. Spectral analysis 

VHE 7-ray spectra result from particle acceleration processes 
and thus they are expected to steepen above a given energy; 
this combines with the energy resolution currently achieved 
by imaging Cherenkov atmospheric detectors (20% at best) to 
cause a considerable event flow into higher estimated energy 
intervals. Starting with an observed differential 7-ray trigger 
rate, one therefore needs a global forward-folding method, us- 
ing the knowledge of the detector response (7-ray effective de- 
tection area, energy resolution), as well as a parameterization 
of the spectral shape. Therefore, we have chosen a maximum 
likelihood method which directly provides relevant physical re- 
sults for the present problem, namely the values of the most 
probable spectral parameters and their covariance matrix. 

The image analysis described in Sect. also yields 
the energy of each hypothesised 7-ray shower. The spectral 
analysis presented below involves the exact energy-resolution 
function T, which is characterised by a rm.s. of 22% (inde- 
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Fig. 3. Mkn 421 nightly-averaged integral flux above 250 GeV between December, 1996, and June, 2000. The 7-ray effective 
area has been weigh ted u sing a differential index of —2.9, in order to estimate the integral flux for observations far from the 
Zenith (see appendix A_^). Arrows stand for 2 a upper-limits when no signal was recorded, and dashed lines show the mean flux 
for each observation year. 



pendent of energy) and which includes possible bias in energy 
reconstruction close to the detection threshold. This function 
has been determined by detailed Monte-Carlo simulations of 
the telescope response, as has the effective detection area AcS, 
which includes the effect of event-selection efficiency (see 
Fig. The simulations have been checked and calibrated on 
the basis of several observables, especially by using muon 
rings and the nearly-pure 7-ray signal from the highest flare of 
Mkn 501 in April 1997 (|Piron et al. 1999a|). 



With typical statistics of '^1000 7-ray events and signal- 
to-background ratio of ^^0.4 (as obtained on the Crab nebu- 
la), a spectrum can be determined with reasonable accuracy 
as follows. First we define a set {A,_ }={[6lf;"\ dT,'''']}^,=un, 
of riz zenith angle bins, with a width (between 0.02 and 
0.04 in cosine) small enough compared to the variation 
scale of T and AeS', Ai corresponds to the transit of the 
source at Themis, and A„, to the maximum angle fixed 
by the data sample. Then we define rig estimated energy 
bins {A,^}={[Ef^''\ Ef^'''']}^^=i^n^, with a width (>0.2 in 
logj^Q -Etcv) 31 least twice as large as the typical width of the 
function T. The maximum energy E™'^^ is fixed by the availa- 
ble statistics. Finally, we define a set of bins {A^^ i^}={Aj^ (g) 
\Ji,=l.n,■A,=^,li^,),n/, for cach A^^ bin, the lowest energy 



(and thus the bin Ai^j(i_)) is determined by the telescope de- 
tection threshold which increases with zenith angle (see Fig. H). 

Within each A^, 2D-bin, the number of events passing 
the selection cuts is determined separately for all ON and 
OFF-source data, and the maximum-likelihood estimation of 
the spectral paramete rs is performed following the procedure 
detailed in appendix A. 1 . The likelihood-function expression 
does not rely on a straightforward "ON-OFF" subtraction as in 
usual spectral analyses, but on the respective Poissonian distri- 
butions of ON and OFF events. In particular, this allows pos- 
sible low statistics to be treated in a rigorous manner No hy- 
pothesis is required on the background (OFF) shape, but two 
hypotheses are successively considered for the differential 7- 

ray spectrum — ^: /) a simple power law, 0q' E^Jy (hyp. Ti-^^), 
dE 

which is often a good approximation, at least within a restricted 
energy range (over one or two orders of magnitude), and //) a 
curved shape, ^-h'^+P"" log^o Et.v) ^^^^ j^csy j^ie lat- 
ter parameterization, previously used by the Whipple group for 
the study of Mkn 421 and Mkn 501 ( jKi-ennrich et al. 1999a| ), 
corresponds to a parabolic law in a log(j/F(t/)) vs. log(i/) re- 
presentation, where i'F{v) = E^—^ and E = hv. 

uE 

The relevance of T^p' with respect to TC^^ is estimated from 
the likelihood ratio of the two hypotheses, which is defined as 



A = 21og 
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one degree of freedom and permits the search for possible spec- 
tral curvature. For each data sample, the spectral law finally 
retained is given by the most relevant parameterization of the 
differential spectrum. In the following, we chose to represent 
each spectrum as a function of the true photon energy by an 
area corresponding to the 68% confidence level contour given 
by the likelihood method. 

3. Results 

3. 1. Data sample and light curves 

The complete data sample consists of observations taken be- 
tween December, 1996, and June, 2000. During these periods, 
the source was systematically observed in a range of zenith an- 
gle extending from close to the Zenith up to 45°. The inten- 
sity of the source did not influence the observation strategy. 
However, a selection based on criteria requiring clear moon- 
less nights and stable detector operation has been applied: this 
leaves a total of 139 hours of on-source (ON) data, together 
with 57 hours on control (OFF) regions. The different light 
curves of the four observation periods are shown in Fig. |^. We 
used a differential index of —2.9, which is repr esentative of all 
spectral measurements presented in Sect. 3^2, to estimate the 
integral flux above 250 GeV for all data, especially those taken 
far from the Zenith: this procedure is detailed in appendix A. 2 . 

As can be seen in Fig. ^ the flux of Mkn 421 changed 
significantly between 1996-97 and 1997-98: almost quiet 
during the first period (with a mean flux $>250 GcV=2.5 ± 
0.9 X 10~^^ cm~^ s~^), the source showed a higher mean 
activity during the second period (<i>>250 GcV=6.0 ± 0.6 x 
10~^^ cm~^ s~^), with small bursts in January and March 
sometimes showing up in excess of the steady flux from the 
Crab nebula (which i s $'^^=14 10 ± 0.35 x 10"" cm^^ g-i 
above 250 GcV, see fPiron 200C| ). In 1998-99, the mean VHE 
emission of Mkn 421 (<I>>250 gcV=2.7±0.9 x lO-^cm-^s-^) 
decreased to a level comparable to that of 1996-97. In spite 
of some activity detected during the winter, the weather con- 
ditions in Themis caused a very sparse source coverage. 
Nevertheless in the beginning of 2000 Mkn 421 showed a re- 
markable increase in activity, exhibiting a series of huge bursts. 
As seen in Fig. ^, the bursts recorded in January and February 
2000 clearly appear as the highest ever seen by CAT from this 
source in four years with a nightly-averaged integral flux cul- 
minating at ~2 ^™ and a large night-to-night variability. 

VHE intra-night variability was also observed on a few oc- 
casions. For instance during the night of January 12-13, the 
source intensity increased by a factor of 3.8 in ^2 hours, from 
~0.6 ^™ to -2.3 (with a per d.o.f of 2.5 for the 
absence of any variation), as can be seen in Fig. ^ (upper-left 
panel). At the bottom of this figure, Mkn 421 light curves are 
also shown for three nights from the 3"^ to the 5**^ February. 
While the fluxes recorded by CAT during the first and last 
nights were stable, respectively $>250 GcV— 1-3 (over 
4 hours) and $>25o GoV— 0-7 (over 1 hour), the source 
activity changed dramatically in a few hours during the second 
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night (February 4-5). The CAT telescope started observation 
while the source emission was at a level of 5.5 This 
flux is comparable to the historically highest TeV flux ever 
recorded, i.e., that of Mkn 501 during the night of April 16*'', 
1997 (Djannati-Atai et al. 1999). In spite of the low source ele- 
vation (6*2=44°) 124 7-ray events with a signal significance of 
6.5 (T were detected during the first 30 minutes of observation. 
This may be compared to the 838 7-ray events and significance 
of 13.8cr obtained during the whole night (see Fig. |l|). After this 
first episode, the source intensity was reduced by a factor of 2 
in 1 hour and by a factor of 5.5 in 3 hours. In Fig. ^, each point 
stands for a —30 min observation but a finer binning in time 
does not show any additional interesting features, confirming 
that CAT started observation after the flare maximum. 



3.2. 1998 and 2000 time-averaged spectra 

The data used in this section consist of a series of —30 min 
acquisitions for which a 7-ray signal with significance greater 
than 3(7 was recorded, and they have been further limited to 
zenith angles 0z<28°, i.e., to a configuration for which the de- 
tector calibration has been fully completed. The spectral study 
is thus based on 6.2 hours of on-source (ON) data taken in 
1998 and 8.4 hours in 2000. Though this data selection reduces 
somewhat the total number of 7-ray events, it provides a high 
signal-to-noise ratio, minimizes systematic effects, and allows 
a robust spectral determination. Concerning systematic effects, 
another favourable factor is the low night-sky background in 
the field of view due to the lack of bright stars around the 
source. 

Systematic errors are thus mainly due to the uncertainty 
on the absolute energy scale, which comes from possible 
variations of the atmosphere transparency and light-collection 
efficiencies during the observation periods. To a lesser extent, 
they are also due to limited Monte-Carlo statistics in the 
determination of the effective detection area. These errors, 
assumed to be the same for all spectra, are implicitly con- 
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Fig. 5. Mkn 421 time-averaged spectra between 0.3 and 5.0 TeV in 1998 and 2000, for the power-law and curved shape hy- 
potheses. The areas show the 68% confidence level contour given by the likelihood method. In each of the four panels, the two 
lower plots give the ratio, in each bin of estimated energy, of the predicted number of events to that which is observed both for 
the 7-ray signal {Ity) as well as for the hadronic background (i?h)- 
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Table 1. Characteristics of the Mkn 421 spectra obtained in this paper. For each spectrum we indicate the observation period, 
the total energy band used in the likelihood method, AE^, the total observed number of 7-ray events, S^, the spectral parameters 
obtained in the T^p' and T-C^^ hypotheses, and the likelihood ratio A. We also quote the decorrelation energy in the T^p' hypothesis, 
Ed, the energy Eq^ at which the energy-dependent exponent 'yf^{ETcv)=l'^^+f3'^^ logio ^TcV has a minimal error in the T-C^^ 



hypothesis (see appendix B.2), and the corresponding value Jq'^ = "ff'^{EQ^). The energies are given in TeV, and the flux 
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7=" 


Ed 


■i>0 




,3" 


i5- 


-'0' 


A 


1998 


0.3 


-5.0 


735±57 


2.29±0.20 


2.88±0.12 


0.69 


2.51±0.32 


2.94±0.15 


0.40±0.48 


0.73 


2.88±0.14 


0.75 


2000 


0.3 


-5.0 


1424±71 


2.90±0.18 


2.95±0.08 


0.63 


3.34±0.28 


3.10±0.13 


0.74±0.37 


0.61 


2.94±0.10 


4.90 


4-5/02/2000 


0.3 


-2.0 


609 ±40 


3.16±0.34 


2.82±0.15 


0.57 


3.25±0.37 


3.08±0.39 


0.63±0.81 


0.37 


2.81±0.17 


0.61 



sidered in the following and they have been estimated from 
detailed simulations ( ^iron 2000| ): (A(/)o/(?!)o)'^''=±20%, 
(A7)"y'^'=±0.06, and (A/3)"y"*=±0.03. 

The 1998 and 2000 time-averaged spectra are shown in 
Fig. ||, both in the power-law and curved shape hypotheses. The 



statistics used for their extraction are detailed in appendix B.l 
the spectral parameters are summarized in Table [l|, and their 
covariance matrices are given in appendix |B.2[ In each panel 



of Fig. H, the two lower plots give the ratio, in each bin of esti- 
mated energy, of the predicted number of events to that which 
is observed both for the 7-ray signal (i?-y) as well as for the 
hadronic background (i?h)- This is another means to check the 
vaUdity of the parameters estimation, and to compare between 
the two hypotheses on the spectral shape. 

As can be seen in Fig. ^a, the power law accounts very 
well for the 1998 time-averaged spectrum. The likelihood ratio 
value is low (A=0.75, corresponding to a chance probability 
of 0.39), and the curvature term is compatible with zero 
(^':«=0.40±0.48"*^*). Thus, we find the following differential 
spectrum: 



dE 



(2.29 ± 0.20"*^' ± 0.46"''"*)10 



-11 —2 

cm s 



1 TeV" 



XEr 



-2.88±0.12="^"±0.06''' 
TcV 



On the contrary. Fig. |]d shows some evidence for a 
curvature in the 2000 time-averaged spectrum, with a rela- 
tively high likelihood ratio value (A=4.90, corresponding to 
a chance probability of 0.027). The lower plots of Fig. 
when compared to those of Fig. |^c, directly confirm that 
a curved spectrum provides a better fit to the data. This 
curvature is equivalent to an energy-dependent exponent 
7f''(£;Tcv)==7''+/3''logio^TcV. Table |l| shows the energy 
Eq^ at which the error on 7f''(£'Tcv) is minimal, as well as 
the corresponding value -fQ^ = 7f^(£'o'^); the latter value is in 
all cases very close to that of 7P' (in the absence of curvature), 
showing the consistency of the spectral analysis. Thus, we fi- 
nally retain the following parameterization for the 2000 period: 



-j^ = (3.34 ± 0.28"*^* ± 0.67"^"*) 10 
dE 



xE. 



"3.10±0.13=*'"±0.06=y='-(0.74±0.37=*'"±0.03=y=') log^o Bt 
TeV 



1997/98 



> 
H 

— 2.8 - 
X 




i HEGRA 

• CAT (this work) 

■ Whipple 



2000 



04/02/00 



21 /04/98 



06/95 



07/05/96 



15/05/96 
Ik — 
28/04/00 
(preliminary) 



Diff. flux % (10 "cm \"'TeV ') at ITeV 



Fig. 6. Compilation of s pectral measurements on Mkn 421 
obtained by the HEGRA ( [Aharonian et al. 1999| [Horns 200 1|) 
CAT (this work) and Whipple ([Zweerink et al. 1997 



Krennrich et al. 1999a , 1999b ) experiments between 1995 
and 2000. The results are given in the {0g',7P^} plane of 
spectral parameters in the power-law hypothesis. While the 
dashed lines around CAT points show only the 68% confidence 
level contour (statistical error), the error bars take account, for 
all measurements, of both statistical and systematic errors, the 
latter being conservative if unknown: (Acj)^^ /(t)Q^Y^^*^=±30% 
and (A7P')'^y'''=±0.10 were used where these values are not 
quoted directly. The per d.o.f corresponding to the absence 
of any dependence of the spectral index with flux is 3.5. 



4. Discussion 

4. 1. Comparison with ottier VHE detections 



The energy spectrum of Mkn 421 has been measured 
i n the same energy range by the Whipple Obs ervatory 
dZweerink et al. 1997| ; |Ki-ennrich et al. 19994 |l999b|) and b y 
the HEGRA experiment ( [Aharonian et al. 1999^ [Horns 200l| ). 
All these experiments have fitted a power law to the spectrum 
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and given the values of and of 7P' in several observation 
periods coiTesponding to different levels of activity of the 
source. Our results for the periods 1998 and 2000 are plotted 
in Fig. H together with their results. This figure includes 
both statistical and systematic errors for all measurements. 
Although CAT data tend to indicate some spectral curvature 
in 2000, we used the values of 7P' for this comparison, given 
that this parameter is in any case clearly indicative of the 



Photon energy E (TeV) 
10 



energy dependence of the spectrum (70*^ 



7 



P', see end of 



Sect. 3.2). In this figure we also added the result obtained 
during the single night between 4 and 5 Febr uary, 2000, for 
which the high source intensity (see Sect. 3A_) allowed us to 
extract a spectrum. Since the corresponding errors are larger 
(see Table |l|), this single result is fully compatible with the 
other two. 

Some caution is necessary when comparing results from 
experiments with different systematic errors. As an exam- 
ple, in the case of a curved spectrum, the value of the ener- 
gy Eq'^ (see Table |l]) corresponding to the minimal error on 
the spectral exponent 7^"^ is not necessarily the same for ex- 
periments with different energy thresholds. Possible systema- 
tic effects could be responsible for the dispersion of results 
shown in Fig. ||, since, if one excepts a recent preliminary 
result from HEGRA ( [Horns 200 1| ), HEGRA and CAT data 
favour higher values of 7P', whereas data from the Whipple 
Observatory yield a rather low spectral index. However, expe- 
riments showing the largest difference in spectral indices in the 
figure (Whipple Observatory and HEGRA) are in good agree 



ment on the spectrum of the Crab nebula ( [Hillas et al. 1998 



Aharonian et al. 2000| ). Furthermore, some evidence for a spec- 
tral variability has been already suggested on the basi s of the 
Whipple 1995-96 low-flux data ( [Zweerinket al. 1997| ), with a 
spectral index 7p'=2.96±0.22, but no flux value is quoted in 
the latter reference, precluding any comparison with the results 
discussed here. 

Data from all experiments with the source in both low and 
high states would be desirable in order to exclude that experi- 
mental effects are responsible for the dispersion seen in Fig. ^. 
Besides, as different flares may be intrinsically different from 
each other, the spectral index-flux correlation may not provide 
us with complete information, thus stressing again the need for 
more numerous observations. 



4.2. Comparison with results on Mkn 501 

The second TeV blazar, Mkn 501, underwent a series of 
intense flares in 1997 observed by all Cherenkov telescopes 
in the Northern Hemisphere (Djannati-Atai et al. 1999 and 
references therein). Fig. ^ shows two spectral energy distri- 
butions (SEDs) of Mkn 501, as measured by CAT. The first 
one is the 1997 time-averaged spectrum and the second one 
corresponds to the sing le night of April 16*'\ 1997, during 
which the source reached its highest intensity. The spectral 
variability observed by C AT on the basis of a hardness ratio 
( Djannati-Atai et al. 1999 ) is illustrated here by the shift in 
the peak energy of the SED: 0.56±0.10 TeV for the average 



> 

> 10 



Mkn 501^ 



^ 



Mkn 421 



> 

a) 

-.0 b 

10 H 
-o 

-o 



V (Hz) 

Fig. 7. Comparison of the Mkn 421 spectral energy dis- 
tributions derived in this paper (unfilled areas) with that 
obtained by CAT for Mkn 501 in 1997 (hatched areas) 
dDjannati-Atai et al. 1999| ; ^iron 2000| ). All areas show the 
68% confidence level contour given by the likelihood method 
in the curved shape hypothesis. 



spectrum and 0.94±0.16 TeV for the highest flare. Fig. ^also 
shows the SEDs of Mkn 421 for the two observation periods 
1998 and 2000: both are obta ined with the assumption of a 
curved spectrum (see Sect. 3.2). 



In the framework of the phenomenological unification 
scheme of blazars established by Fossati et al. (1998), the 
broad-band SED of such a source consists of two components, 
respectively a low-energy component attributed to synchrotron 
radiation and a high-energy component peaking in the 7-ray 
range. Peak energies of the preceding components are corre- 
lated and Mkn 501, as observed in 1997, appears to be the most 
extreme blazar, with the synchrotron part of the SED spectrum 
peaking in the hard X-ray range and the 7-ray component pea- 
king around or above 500 GeV. 

In the case of Mkn 421, the CAT 1998 and 2000 data, and 
the weakness of its detection by EGRET, when interpreted 
within the preceding unification scheme, indicate that the 
peak energy of the high-energy component of the SED lies 
in the 50 GeV region, while its synchrotron peak is known 
to lie in the UV-soft X-ray band. A shift of both components 
towards higher energies during an intense flare, similar to that 
observed by CAT for Mkn 501 in 1997, would result in a 
harder spectrum in the energy region covered by Cherenkov 
telescopes and thus account for the trend shown in Fig. ^ 
In fact, the spectral hardening of Mkn 421 at TeV energies 
during flares, if clearly proven, should not be a surprise, 
since a similar effect has been already observed at a few keV 
from this source by t he ASCA and Be ppo-SAX satellites 
( fTakahashi et al. 1999| ; ^alizia et al. 200C| ), and since X-ray 
and 7-ray integrated emissions of Mkn 42 1 have proven to be 
correlated on many occasions (see, e.g., Maraschi et al. 1999; 
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Takahashiet al. 1999, 2000). 



SSC models have been already successfully applied to 
the Be ppo-SAX and CAT data obtained on Mkn 501 in 1997 
(e.g., ^iron 206o| ; Katarzynski et al. 2001 ), but alternative 



5. Conclusions 

Along with Mkn 501, Mkn 421 is the most extreme known 
BL Lac object. The CAT results presented in this paper con- 
firm the high variability of this source at TeV energies, which 
has been already reported in the past for intense and rapid 
bursts (paidos et al. 19961). Whereas Mkn 501 TeV variabili- 



ty has been observed by CAT from night to night, but never 
on shorter time-scales, Mkn 421 showed a fast variation in in- 
tensity within one hour during the night from 4 to 5 February, 
2000. A simple causality argument implies that the 7-ray emit- 
ting region must be very compact with a size of ~10 light-hours 
if one takes a typical value of 10 for the geometric Doppler fac- 
tor (see, e.g., Celotti et al. 1998 ), which reduces the time-scale 
in the observer frame. 

As for the spectral properties, Mkn 421 has proved to 
be a little less extreme than Mkn 501. Whereas the latter 
has shown a 7-ray peak lying above the CAT threshold, 
Mkn 421 exhibited a power-law spectrum in 1998, indicating 
that current imaging Cherenkov detectors cover the end 
part of its spectral energy distribution. However, the 2000 
time-averaged spectrum shows some indication of curvature, 
which in fact has been also marginally observed by the 
Whipple Observatory at the time of the 1995-96 flaring 
periods (Krennrich et al. 1999a). 



The present observations of Mkn 421 leave some open 
questions. Firstly, definite conclusions on a possible spec- 
tral variability still need simultaneous observations of more 
numerous flares by different Cherenkov telescopes in order 
to exclude possible systematic effects between experiments. 
Secondly, the interpretation of the multi-wavelength spectra 
of Mkn 421 still suffers from the lack of precise simultane- 
ous spectral measurements of the low and high-energy parts 
of its SED which could bring accurate constraints to mo- 
dels. For instance, in the Synchrotron Self Compton (SSC) 
model wh ich is often invoked t o explain the SED of extreme 
BL Lac's ( Ghisellini et al. 1998 ), the 7-ray component is inter- 
preted as the result of the inverse Compton process occurring 
between the ultra-relativistic electrons, which emit synchrotron 
radiation at low energies, and the latter soft photon field it- 
self. This simple and most natural model predicts a strong 
correlation between the synchrotron and the 7-ray bump be- 
haviours. Such a correlation has been observed many times on 
Mkn 421, for example during the coordinated observation cam- 
paign in spring 1998, which involved ground-based Cherenkov 
imaging telescopes (Whipple, HEGRA, and CAT), and the 
Beppo-SAX and ASCA X-r ay satelHtes ( [Maraschi et al. 1999 ; 
Takahashi et al. 1999, 200C). The correlation was only proven 
in terms of integrated (and not differential) fluxes due to the 
lack of statistics. Similarly in 2000, the TeV flaring behaviour 
of Mkn 421 was accompanied by an overall increase of its 
keV activity, as for instance continuously recorded by the ins- 
truments on board the Rossi X-Ray Timing Explorer satellite 
dPiron 2000|). 



scenarios still exist also for this source (e.g., Rachen 1999 
Muecke et al. 2001), which consider an e^p plasma with 



synchrotron photons radiated by electrons at low energies and 
7-rays emitted at high energies by the products of the proton- 
induced cascades (from tt'^ decays and from synchrotron 
radiation of protons and muons). Thus, in order to address 
more deeply the crucial problem of the plasma jet content, a 
richer sample of flares from both blazars, detected at various 
wavelengths, would be necessary. 



Acknowledgements. The authors wish to thank the French na- 
tional institutions IN2P3/CNRS and DAPNIA/DSM/CEA for sup- 
porting and funding the CAT project. The CAT telescope was also 
partly funded by the Languedoc-Roussillon region and the Ecole 
Polytechnique. The authors also wish to thank Electricite de France 
for making available to them equipment at the former solar plant 
"Themis" and allowing the building of the new telescope and its 
hangar. They are grateful to the French and Czech ministries of 
Foreign Affairs for providing grants for physicists' travel and accom- 
modation expenses. L.R. thanks for the financial support granted by 
the Ministry of Education of the Czech Republic (Project LN00A006). 

References 

Aharonian, F.A., Akhperjanian, A.G., Andronache, M., et al. 1999, 
A&A 350, 757 

Aharonian, F.A., Akhperjanian, A.G., Barrio, J.A., et al. 2000, ApJ 
539, 317 

Barrau, A., Bazer-Bachi, R., Beyer, E., et al. 1998, NIM A416, 278 
Buckley, J.H., Akerlof, C.W., Biller, S.D., et al. 1996, ApJ 472, L9 
Celotti, A., Fabian, A.C., & Rees, M.J. 1998, MNRAS 293, 239 
Chariot, P, Gabuzda, D., Sol, H., Degrange, B., & Pare, E. 1998, ig**" 

Texas Symposium on Relativistic Astrophysics and Cosmology 

(Paris) 

Djannati-Atai, A., Piron, F, Barrau, A., et al. 1999, A&A 350, 17 
Fossati, G., Maraschi, L., Celotti, A., Comastri, A., & Ghisellini, G. 

1998, MNRAS 299, 433 
Gaidos, J.A., Akerlof, C.W., Biller, S.D., et al. 1996, Nature 383, 319 
Ghisellini, G., Celotti, A., Fossati, G., Maraschi, L., & Comastri, A. 

1998, MNRAS 301,451 
Hartman, R.C., Bertsch, D.L., Bloom, S.D., et al. 1999, ApJSS 123, 

79 

Hillas, A.M., Akerlof, C.W., Biller, S.D., et al. 1998, ApJ 503, 744 
Holder, J., for the CELESTE collaboration 2001, High-Energy 

Gamma-Ray Astronomy (Heidelberg), AIP Conf. Proc. 558, ed. 

FA. Aharonian & H.J. Voelk, 635 
Horns, D., for the HEGRA collaboration 2001, 36**" Rencontres de 

Morio nd on "Very High- Energy Phenomena in the Universe" (Les 

Arcs), |astro-ph/0103514| 
Katarzynski, K., Sol, H., & Kus, A. 2001, A&A 367, 809 
Krennrich, F, Akerlof, C.W., Buckley, J.H., et al. 1997, ApJ 481, 758 
Ki-ennrich, F, Biller, S.D., Bond, I.H., et al. 1999a, ApJ 511, 149 
Ki-ennrich, F, Bond, I.H., Bradbury, S.M., et al. 1999b, 26* 

International Cosmic-Ray Conference (Salt-Lake City) 3, 301 
Le Bohec, S., Degrange, B., Punch, M., et al. 1998, NIM A416, 425 
Lin, Y.C., Bertsch, D.L., Chiang, J., et al. 1992, ApJ 401, L61 
Macomb, D.J., Akerlof, C.W., Aller, H.D., et al. 1995, ApJ 449, 99 
Malizia, A., Capalbi, M., Fiore, F, et al. 2000, MNRAS 312, 123 189 



10 



F. Piron et al.: gamma-ray properties of Mkn 421 above 250 GeV from CAT observations between 1996 and 2000 



Maraschi, L., Fossati, G., Tavecchio, F, et al. 1999, ApJ 526, L81 
Muecke, A., & Protheroe, RJ. 2001, APh 15, 121 
Petry, D., Bradbury, S.M., Konopelko, A., et al. 1996, A&A 311, L13 
Piron, F, for the CAT collaboration 1998, 19"^ Texas Symposium on 

Relativistic Astrophysics and Cosmology (Paris) 
Piron, F, for the CAT collaboration 1999a, GeV-TeV Gamma-Ray 
Astrophysics Workshop (Snowbird), AIP Conf. Proc. 515, ed. 
B.L. Dingus, M.H. Salamon & D.B. Kieda, 308 
Piron, F., for the CAT collaboration 1999b, 26"^ International Cosmic- 
Ray Conference (Salt-Lake City) 3, 326 
Piron, F. 2000, Ph.D. Thesis, University of Orsay-Paris XI 



( tittp://lpnhp2. in2p3.fr/homecat/publi.htm: ) 
Punch, M., Akerlof, C.W., Cawley, M.F, et al. 1992, Nature 358, 477 
Rachen, J.P. 1999, GeV-TeV Gamma-Ray Astrophysics Workshop 

(Snowbird), AIP Conf. Proc. 515, ed. B.L. Dingus, M.H. Salamon 

& D.B. Kieda, 41 
Takahashi, T, Tashiro, M., Madejski, G., et al. 1996, ApJ 470, L89 
Takahashi, T, Madejski, G., & Kubo, H. 1999, APh 11, 177 
Takahashi, T, Kataoka, J., Madejski, G., et al. 2000, ApJ 542, L105 
von Montigny, C, Bertsch, D.L., Chiang, J., et al. 1995, ApJ 440, 525 
Zweerink, J.A., Akerlof, C.W., Biller, S.D., et al. 1997, ApJ 490, L141 

Appendix A: Spectra and light curve extraction 

A.1. The spectral analysis procedure 

Here we give the expression of the HkeHhood function C used 
for s pect ral reconstruction, using the same notations as in 
Sect. 2.3. We start with an assumption on the spectral shape 



by choosing a given law ( — ^ 

dE 



prcd 



which includes the set 



{A} of parameters to fit. Then, for each 2D-bin ^i^^i^, cor- 
responding to the zenith angle interval [6*""", Ofl'^'^] and to the 
estimated energy interval [i?'"'", Ef^^^], let us define: 

- and the numbers of events passing the selec- 
tion cuts in the ON and OFF data, respectively; 

- as the normalisation factor between ON and OFF ob- 
servations, obtained from the contents of the respective a 
distributions in the range from 20° to 120° (the control re- 
gion) (see Fig. |l|); 

~ = "ij.ic ~ PizPi^,ic as the obscrvcd number of 7-ray 

events, and SSi^^i^ 

-fpred 



+ 0f^Pi,.i^ as its error; 



as the predicted number of 7-ray events: if we note 



T(6'z, E^—^Ej) dE~f for the probability, at fixed zenith an- 
gle 9z and real energy E^, to get an estimated energy E^ 
within the interval [E^, E^+dE-^], then 



dE 



dE 



dE 



prcd 



AM^,E)r{e,^,E-^E), 



(A.l) 



where Tqn is the total ON-source time of observation, ^eff 
is the effective detection area, and where is defined by 
cos(5~)=i [cos(0f"") + 008(6'^;;^'^)] ; 

- Pi~i^ the mean predicted number of events in the OFF data; 

- tmTJI = Sf^'f + Pi^PiTJI the mean predicted number of 
events in the ON data. 



The observed numbers and Pi^^i^ have Poissonian 

probability distributions V{ni^^i^) and 'P{Pi^,i^), respectively, 
and the likelihood function is as follows: 

^{{^},{w:j:}) - n ^("«..^e)^(K..^J- (A.2) 



The quantities Pi,.i^, which are unknown, can be 

determined by maximizing the function C: the re- 

d\og(C) 

lations =0 lead to the solutions Pi^,i^ = 



1 



2/3^. (A. +1) 
with fli^^j^ = PiA'^i., 



prcd 



prcd 



Finally, we reinject these expressions in Eq. ( |A.2| ), elimi- 
nate the "constant" terms which only depend on or pi^^i^, 
and get: 



log [£ ({A})] = h^-^= i°g (Ct + 



S: 



prcd 



(A.3) 



The { A} parameters, included in this expression through 
Eq. ( |A.1| ), are the only unknown quantities left: their values 
are determined using an iterative procedure of maximization, 
which leads to the final fitted values and their covariance matrix 
V. 

A.2. Integral flux determination 

When computing an integral flux, one must again take care of 
the telescope detection threshold increase with zenith angle O^; 
thus, let us denote: 

- £""'[6*2] for the unique energy which verifies 

Aff (0z, ^'"*) = Aff(0°, 250 GeV); 

- S*'"* for the total number of 7-ray events observed within 
the selection cuts with an estimated energy above E™^[9^ 
(S*'"* is determined by a simple "ON-OFF" subtraction), 
and (55*'"* for its error 

To convert S*'"' into the corresponding integral flux $ above 
250 GeV, we use the following proportionality relation: 



^iiit, best 
best 



250 GcV 



best 



where ( ) is the source differential spectrum measured 



as explained in Sect. 23, and S*'"*- is the number of 7-ray 
events predicted from this law above i?'"*[0z]: 



S 



int. best 



ON 



dE 



dE 



dE 



best 



A,s{0,,E)T{e,,E-^E). 



Finally, the error on $ is simply estimated as (5$ — —^r^ x $. 
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13.1 
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300 261 
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10.4 
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1.3 
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26 


5.7 
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411 
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30 


9.2 


300 


172 128 


18 


7.0 


1.3- 


2.0 


229 


123 


107 


18 


5.9 


335 


158 


175 


21 


8.6 


138 
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12 


6.0 


2.0- 


5.0 


119 


66 


53 


13 


4.0 


165 


89 


76 


15 


5.1 











Table B.2. Covariance matrices elements of Mkn 42 1 spectra. 



Period vg v|| 

1998 3.97x10-2 -1.24x10^2 1.44x10"^ 1.03x10-^ -1.10x10"^ 1.13x10-^ 2.35x10-^ 3.28x10"^ 2.35x10" 

2000 3.19x10^2 -8.65x10"'' 6.51x10"-' 7.80x10"^ -8.36x10"* 6.45x10"^ 1.68x10"^ 2.91x10"^ 1.36x10" 

4-5/02/2000 1.12x10"! -3.84x10"^ 2.16x10"^ 1.38x10"^ -3.45x10"^ 3.63x10"^ 1.53x10"^ 2.86x10"^ 6.57x10" 



Appendix B: More on Mkn 421 spectra 

B. 1. Available statistics 

The statistics used for extracting the Mkn 421 spectra are given 
in Table [B.l[ With the same notations used in this paper, we 
write: 

- the data zenith angle amplitude, and the ON and OFF- 
source observation durations. Ton and Toff; 

- for each of the rig energy intervals A^^ = [i^™'", Ef^'^^] (in 
TeV): 

• the numbers of events passing the selection cuts in the 
ON and OFF data, 

• the observed number of 7 events, Si^ = n^^ — pi^ , and 

its eiTor, SSi^ = {Y,^^ rii^^i^ + Pi.'^Pr.sS ^' 

• the statistical significance N^i^ — Si^/SSi^; 

- the observed total number of 7 events, = J^i ' ^^'^ 

its error, ,55-, = {J2, SSi^^)^^^. 



At this decorrelation energy, the width of the hatched area re- 
presenting the 68% confidence level contour, is minimum: for 
instance, we find £'d=690 GeV for the 1998 time-averaged 



spectrum of Mkn 421 shown in Sect. 3.2 



In the curved shape hypothesis, the energy-dependent ex- 



ponent 7f" (E'toV ) 



at the energy E^"" 
•jf^ writes 



^cs _|_ ^cs log-^Q -EtoV has a minimal error 
10 ^f^f TeV; the corresponding value of 



7- ^ jf^iE^n = 7- - 



7/3 



V' 



1313 



and its eiTor 



v: 



7/3 



(313 



B.2. Covariance matrices 

Let and V^^ be the two covariance matrices o f sp ectral 
parameters for the TC^^ and H'^^ hypotheses (see Sect. 2.3 ), res- 
pectively: 



/ T/Pl T/Pl \ 

ypi ^ 4'<l> <t>i and 

\ 07 77 / 



T/CS T/CS T/CS 

"^4,4, ^<j,f3 

T/CS T/CS T/CS 

*^07 ^77 *^7/3 

T/CS T/CS T/CS 

^4,0 ^^13 ^fjp . 



The values of these matrices elements are given in Table B.2 
for each spectrum of Mkn 42 1 obtained in this paper 



In the power-law hypothesis, the energy E^ at which the 
values of (hf} et 7P' are decorrelated writes: 



Ea 



exp 



pi 



4>f 



xplypl 



TeV. 



